Distribution of chirality in the quantum walk: Markov process and entanglement 
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The asymptotic behavior of the quantum walk on the line is investigated focusing on the probabil- 
ity distribution of chirality independently of position. It is shown analytically that this distribution 
has a long-time limit that is stationary and depends on the initial conditions. This result is un- 
expected in the context of the unitary evolution of the quantum walk, as it is usually linked to a 
Markovian process. The asymptotic value of the entanglement between the coin and the position is 
determined by the chirality distribution. For given asymptotic values of both the entanglement and 
the chirality distribution it is possible to find the corresponding initial conditions within a particular 
class of spatially extended Gaussian distributions. 
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I. INTRODUCTION 

The quantum walk (QW) on the line [l| is a natural 
generalization of the classical random walk in the frame 
of quantum computation and quantum informatkin pro- 
cessing and it is receiving permanent attention [2- 4]. It 
has the property to spread over the line linearly in time as 
characterized by the standard deviation a{t) ^ t, while 
its classical analog spreads out as a{t) ~ t^/^. This prop- 
erty, as well as quantum parallelism and quantum entan- 
glement, could be used to increase the efhciency of quan- 
tum algorithms. As an example the QW has been used as 
the basis for optimal quantum search algorithms [5, 6] on 
several topologies. On the other hand, some experimen- 
tal implementations of the QW have been reported 0, 
and others have been proposed by a number of authors 

i- 

The concept of entanglement is an important element 
in the development of quantum communication, quan- 
tum cryptography and quantum computation. In this 
context, several authors have studied the QW subjected 
to different types of coin operators and/or sources of de- 
coherence to analyze the long-time entanglement between 
the coin and the position and its relation with the ini- 
tial conditions. Carneiro et al. investigate entangle- 
ment between the coin and the position calculating the 
entropy of the reduced density matrix of the coin. The 
relation between asymptotic entanglement and nonlocal 
initial conditions (in the one and two dimensional QW) 
is treated in p^| - [l7| . Ref. [T8l - [20j analyze the effect of 
entanglement on the initial coin state, which is measured 
by the mean value of the walk, and the relation between 
the entanglement and the symmetry of the probability 
distribution. In Ref. [2l| the relation between entan- 
glement and decoherence is studied numerically. Refs. 
[2^ [23| propose to use the QW as a tool for quantum 



algorithm development and as entanglement generator, 
potentially useful to test quantum hardware. 

In the previous works [ll| - |23j the QW evolution was 
studied using the amplitude of probability to evaluate 
the dynamics. In this work we introduce a new probabil- 
ity distribution, the global chirality distribution (GCD) 
that is the distribution of chirality independently of the 
walker's position. We show that the GCD has an asymp- 
totic limit and we connect this limit with the entropy 
of entanglement between the coin and the position. The 
asymptotic behavior of GCD is an unexpected result be- 
cause, due to the unitary evolution, the QW does not 
converge to any stationary state as would be the case 
e.g. of a Markov chain. In order to show these results, 
we rewrite the QW evolution equation as the sum of two 
different terms, one responsible for the classical diffusion 
and the other for the quantum coherence 2J]. As we 
shall see, the first term above obeys a master equation 
as is typical of Markovian processes, while the second 
term includes the interference needed to preserve the uni- 
tary character of the quantum evolution. This approach 
provides a more intuitive framework which proves use- 
ful for analyzing the behavior of quantum systems with 
decoherence. It allows to study the quantum evolution 
together with the associated classical Markovian process 
at all times, and in particular the asymptotic behavior of 
the GCD. 

The paper is organized as follows. In the next section 
we develop the standard QW model, in the third section 
we build the master equation for the GCD, in the fourth 
section we present the asymptotic solution for the QW, 
in the fifth section the entropy of entanglement is con- 
nected with the GCD, and in the last section we draw 
the conclusions. 



II. THE STANDARD QW 
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The QW on the line, corresponds to a one-dimensional 
evolution of a quantum system in a direction which de- 
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pends on an additional degree of freedom, the chirality, 
with two possible states: "left" \L) or "right" \R). The 
global Hilbert space of the system is the tensor product 
Hs ® He where Hs is the Hilbert space associated to the 
motion on the line and He is the chirality Hilbert space. 
Let us call T- (r+) the operators in Hs that move the 
walker one site to the left (right), and \L){L\ and |-R)(-R| 
the chirality projector operators in H^- We consider the 
unitary transformations 



where K{0) = Oz cos9 + iax sin 6*, / is the identity oper- 
ator in Hs, and az and ax are Pauli matrices acting in 
He- The unitary operator U{9) evolves the state in one 
time step as 

\^{t+i)) = umnt))- (2) 

The wave vector can be expressed as the spinor 



k— — oo 



ak{t) 
hk{t) 



\k), 



(3) 



where the upper (lower) component is associated to the 
left (right) chirality. Substituting Eq. ^ and Eq. (P) 
into Eq. ^ and projecting over the position vector \k) 
the unitary evolution is written as the map 

ak{t+l) = ak+i [t] cos 6 + bk+i{t) sin 6*, 
bk{t+l) = ak-i{t) sine - bk^i{t) cos9. (4) 



III. UNITARY EVOLUTION AND MASTER 
EQUATION FOR THE CHIRALITY 

In references (23 . [2H , it is shown how a unitary quan- 
tum mechanical evolution can be separated into Marko- 
vian and interference terms. Here we use this method 
to recognize a master equation in chirality starting from 
the original map Eq. First we define the left and 

right distributions of position as Pkhit) = \0'k{t)\^ and 
PkR{t) = 1 6fc(t) I ^ respectively. Combining the two com- 
ponents of the Eq. ^ and after some simple algebra we 
obtain 

PkAt + 1) = Pk+iAt) cos"" e + Pk+iMt) sin' 
+ I3k+i{t) sin 20, 

PkMt + 1) = Pk-iAt) sin' + Pk-iMt) cos' e 

~ I3k-i{t) sin2e , (5) 

where f3k = Re [ak(t)bl.{t)] is an interference term, with 
Re(z) indicating the real part of z. Of course the prob- 
ability distribution for the position is Pk{t) = Pkhit) + 
PkR{t)- We define the global left and right chirality prob- 



abilities as 

00 00 

PUt)^ E E 

fc— — 00 k— — oo 

00 00 

Pn{t)^ E E I^'^WI'' 



(6) 



k^-c 



k— — c 



Uie)^{T^®\L){L\ + T+®\R){R\}o{I®K{0)}, (1) couple 



with Pii{t) + PL{t) = 1 and the global chirality distribu- 
tion (GCD) is defined as the distribution formed by the 

PL{t) 



PR(.t) 



Using the definition Eq. © in Eq. ([5]) we have 



PUt + i) 

PR{t + 1) 



cos^ 9 SIT? t 
sin^ 9 cos^ I 



Re[Q(t)]sin26l 



PUt) 
Pnit) 



where 



(7) 



(8) 



In Eq. ([7]) the two dimensional matrix can be interpreted 
as a transition probability matrix for a classical two di- 
mensional random walk as it satisfies the necessary re- 
quirements, namely, all its elements are positive and the 
sum over the elements of any column or row is equal to 
one. On the other hand, it is clear that Q{t) accounts for 
the interferences. When Q{t) vanishes the behavior of 
the GCD can be described as a classical Markovian pro- 
cess. However Q{t) — does not necessarily imply the 
loss of unitary evolution; such a loss requires the vanish- 
ing of ah the /3fc(t) (see Eq ([5])). As shown in Ref. ^ 
the primary effect of decoherence is to make the interfer- 
ence terms Pk{t) negligible; in this case Eq. ^ becomes 
a true master equation. On the other hand, when Q{t) 
is time independent, that is Q{t) = Q = constant, then 
Eq. ([7]) is solved using the methods developed in [2^ ; its 
solution as a function of the initial GCD is 



Pdt) 

PB{t) 



1 / 1 + COS* 29 1 - cos* 29 

2 I 1 - cos* 26* 1 -I- cos* 29 



Pl{0) 

Pr(.o) 



+Re [Q] 



1 - cos* 29 
tan 9 



1 

-1 



(9) 



Taking the limit t — > cx) in Eq. (|9]) it is possible to obtain 
the asymptotic value of the GCD as a function of its 
initial value and Q. 

In the generic case Q{t) is a time depend function but 
in this system (as will be seen in the next section) Q{t), 
Phit) and Puit) have long-time limiting values which are 
determined by the initial conditions of Eq. Therefore 
we can solve Eq. ^ in this limit, defining 



Hi = PL{t^ 00), 
Hfl = Pfl(t^oo), 
Qo = Qit 00), 



(10) 
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and substituting these asymptotic values in Eq. ([7]), to have 
obtain the stationary solution for the GCD 



"Hi ■ 


1 




^ 2 



1 + 2Re((3o)/tan( 
1 - 2Re(Qo)/tan( 



(11) 



Qit)= J2 4bt^Qo, 



k— — oo 



(13) 



This interesting result for the QW shows that the long- 
time probability to find the system with left or right chi- 
rality has a limit. 

In the next section we show that it is possible to have 
Qo = choosing adequately the initial conditions. In 
this case, Eq.© approaches a Markov chain Q with 
two states and the dynamics of the GDC turns into an 
example of dependent Bernoulli trials in which the prob- 
abilities of success or failure et each trial depend on the 
outcome of the previous trial. Now the only asymptotic 
solution is Hl = = 1/2 (see Ea. (|TT|) V 

If we look back at Eq. ^ in connection with Eq. ([TT|) a 
paradoxical situation arises. The dynamical evolution of 
the QW is unitary but the evolution of its GCD has an 
asymptotic value characteristic of a diffusive behavior. 
This situation is further surprising if we cornpare our 
case with the case of the QW on finite graphs [23| where 
it is shown that there is no converge to any stationary 
distribution. 



PL(t) 



j=-oo 



k— — oo 



(14) 



(15) 



The time independence of Q{t), PL{t) and PR{t) is a con- 
sequence of the asymptotic approach given by Eq. (jl2p 
and evidently their values are Qo, Hl and II ji respec- 
tively. When Q{t), PL{t) and PR{t) are calculated with 
the map Eq. they have a transient time dependence 
(for t < to) after which they attain their asymptotic val- 
ues Qo, Hi and Hr, as shown in (25j . 

We propose for the initial conditions a° and 6^ the 
following extended Gaussian distributions [29| 



4^ 



1 



(Tov 27r 



exp 



{k~kof 



2al 



cos a, 



(16) 



IV. ASYMPTOTIC SOLUTION FOR THE QW 

In previous works [3, HEl an alternative analytical ap- 
proach was presented to obtain the asymptotic behavior 
of the QW on the line. The discrete map was substi- 
tuted by two continuous differential equations for ak{t) 
and bk{t) starting from a characteristic time » 1. 
The initial conditions for these equations are not neces- 
sarily the same as those used in the discrete map Eqs. 

, because the approximation of a finite difference by a 
derivative does not hold for small times. However these 
initial conditions must assure the same asymptotic be- 
havior than that of the discrete map. 

The asymptotic solutions of Eqs. (U) given by the dif- 
ferential equations are 



ak{t)^ E Jk-i{t cose), 

I — — CO 

oo 

bk{t)^ E (-1)'"'^? Jk-i{t cose), (12) 



where J; is the Ith order cylindrical Bessel function and 
a1 and are initial amplitudes for the differential equa- 
tions. To secure that the behavior of the discrete map 
and the differential equations are the same in the asymp- 
totic regime we should choose a° and b^ to be smoothly 
extended in space. 

Replacing Eq. (IT2|) in Eqs. dHH]) and noting that the 
Bessel functions satisfy X]j^-oo ~ we 



4 



tan a exp(i(5), 



(17) 



where ctq is the initial standard deviation, fco is the central 
position of the Gaussian distribution, a is a parameter 
that determines the initial proportion of the left and right 
chiralities and 5 is a phase to be determined bellow as a 
function of a and 9. 

Now we evaluate Qq, Hl and Hr using the initial con- 
ditions Eqs. (mini) in Eqs. (HI [Til [15]) and noting that 



k=- 



, exp 



27rCTo, for (To > 1 



Qo = 2 ^'^ '"'^^ ^ ' 



cos a , 



siv? a . 



(18) 



(19) 



(20) 



On the other hand, from Eq. (|TT|) we see that Qo and 11^ 
are not independent, then substituting Eqs. (fT8l[T9l[20|) 
into Eq. ([TT|) we have 



cos (5 = 



tan^ 
tan 2a 



(21) 



Then, we rewrite Eq. (IT51) as a function of the two inde- 
pendent parameters of the model e and a 



?o = ^ cos 2a tan0 , 



(22) 
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FIG. 1: The global left probability Phit) as a function of the 
dimensionless time t calculated with the map Eq. Q. The 
initial conditions are given by Eqs. (|16lll7p . It is shown for 
three values of 0, ao — 100, fco = and cos^ a — 0.3. The ap- 
proximate values of to are, from top to bottom, 300, 500, 700. 
The asymptotic values for the dimensionless entropy So are 
also presented. 



note that Qo vanishes for a = 7r/4. In order to ver- 
ify the approximations made in our analytical treatment 
we shall compare the result of Eq. ([T^ with the nu- 
merical evaluation of the asymptotic behavior of Phit) 
using the map Eq. Q and the initial conditions given 
by Eqs. (|16[|17p . This calculation is presented in Fig. [T] 
Our selection of the initial amplitudes makes that the 
asymptotic value for PL{t) is equal to the initial one. 
Which in turn is the same as that given by Eq. (IT^ . 
Thus the asymptotic behaviors of Eq. @ is in excellent 
agreement with our theoretical approach. Our treatment 
works very well for values of ctq > 10. The asymptotic 
regime of Pl {t) sets in at the time after some strong 
oscillations. The value of depends on the parameters 
of the problem. 

To sum up, the dynamical evolution of the QW is de- 
fined by Eq. but it is possible to obtain a prede- 

using as 



termined asymptotic value of the GCD 



initial conditions Eqs. (116117^ . where the parameters a 
and 5 are determined by the Eqs. (|19l2ip and the pa- 
rameters fco and ctq are free to be adjusted. 



V. ENTROPY OF ENTANGLEMENT 

The unitary evolution of the QW generates entangle- 
ment between the coin and position degrees of freedom. 
This entanglement will be characterized [ill by the 
von Neumann entropy of the reduced density operator, 
called entropy of entanglement. The quantum analogue 
of the Shannon entropy is the von Neumann entropy 



Sn{p) = -tr(plogp). 



(23) 



where p = \^ (t)) {"^ {t)\ is the density matrix of the quan- 
tum system. Due to the unitary dynamics of the QW 
the system remains in a pure state and this entropy van- 
ishes. However, for these pure states, the entanglement 
between the chirality and the position can be quantified 
by the associated von Neumann entropy for the reduced 
density operator 



S{p) = -tr(pc log pc). 



(24) 



where pc = tr(p) and the partial trace is taken over the 
positions. Using the wave function Eq. ([3]) and its nor- 
malization properties, the reduced density operator is ex- 
plicitly expressed as 



_( PL{t) Q{t) 
P'-\ Qity Pn{t) 



(25) 



The reduced entropy can be expressed through the two 
eigenvalues {A4., A_} of the reduced density matrix as 



Sit) 



-A+ log2 A+ - A_log2 A_. 



(26) 
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FIG. 3: Level curves for the dimensionless entropy of entan- 
glement So as a function of the dimensionless angle ^ and 
the asymptotic global left probability IIl. Five curves, in full 
line, are presented; each curve has two branches placed sym- 
metrically on both sides of the central straight dashed line 
TIl ~ 0.5. Starting from this line, where So ~ 1, the values 
of So are 0.99, 0.95, 0.90, 0.85 and 0.70. 



The expressions for the eigenvalues are 



l±,/l+4 |(3(t)|'-Pi(()ft(f) 



In the asymptotic regime X± — > A± where 



(27) 



1± Wl + 4 IQol -Hi 



and the corresponding entropy {S(t) — ^ Sq) is 

5*0 = -A+ log2 A+ - A_ log2 A_. 

Using the initial conditions Eqs. p^[T7)) in Eq. 
have 



1 ± 



COS 2a 
cos 9 



(28) 



(29) 



we 



(30) 



For a — tt/A both eigenvalues are A± = 1/2 and from 
Eqs. (Ullinillll) Hi = Hr = 1/2 and Qo = 0. For 
this value the entropy of entanglement Eq. (|29p has its 
maximum value So — 1. Therefore the maximum value of 
the entropy of entanglement is achieved for the classical 
Markovian process (Qo = 0). Note that this result is 
true for all initial conditions that satisfy Q{t) — >■ 0, as it 
follows from Eqs. ©[IIlllHllSni) . 



For a = 9/2, the entropy attains its minimum value 
So = 0, see Eqs. (|291l30p . Then in this case there is no 
entanglement between coin and position. 

Using the results of the previous section it is clear that 
starting from given initial conditions the asymptotic val- 
ues Hl and Qo are obtained, and then the entropy of 
entanglement is calculated using Eqs. (PSIE^ . The in- 
verse path is also possible, that is starting from a prede- 
termined value of the entropy of entanglement Eq. 
it is possible to obtain the initial conditions Eqs. p^flT)) 
of the system that produce this entanglement asymptot- 
ically. 

The previous ideas are numerically implemented using 
Eqs. (ITTIBSIB^ and taking Qo as a real constant, and the 
results are presented in Figs. [2] and [3l Fig. [2] shows that 
for each value of 5*0 there are two values of 11^ and that 
the width of the entropy curve grows inversely with 9. In 
Fig.|3]the level curves of the entropy as a function of 6* and 
Hl are presented as projections of the three-dimensional 
surface. From these figures it is clear that the maximum 
of the entropy of the entanglement is achieved for the 
classical Markovian process (11^ — 1/2) 

To conclude this section, it is interesting to compare 
the entropy of entanglement with the usual Shannon en- 
tropy, used in the theory of communication. In particular 
one could wonder if the entropy of the entanglement may 
be used as a measurement of the degree of disorder of 
chirality. The Shannon entropy, in the asymptotic GCD 
model, is 



Ss = -Ul log2 Ul - Ur log2 Ur, 



(31) 



where 11^ and Ur are given by Eq. ([TT|) . 

It is clear from Eqs. d^HEIl) that when Qo = 
(Hi = Ur = 1/2) both entropies attain the maximum 
value So = Ss = I- However for other values of Qo they 
are different, in particular, when there is perfect statis- 
tical order, i.e. Ul = 1 and Ur = (or 11^ = and 
Ur = 1), the Shannon entropy vanishes as it should but 
the entropy of entanglement does not vanish. Therefore, 
although the behavior of the entropy of entanglement is 
correlated with the behavior of the GCD, it does not de- 
scribe correctly the degree of disorder of the GCD. 



VI. CONCLUSION 

This work provides a different insight for the QW dy- 
namics. It studies the QW focusing on the probability 
distribution of the chirality independently of the posi- 
tion (GCD) and connects this distribution with the en- 
tropy of entanglement. Using an alternative analytical 
approach for the QW on the line, developed in previous 
works [1, [2^, we show analytically that the GCD con- 
verges to a stationary solution. The asymptotic behavior 
of the GCD looks like the behavior of the two dimen- 
sional classical random walk but, unlike the latter, the 
asymptotic GCD depends on the initial conditions. The 
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coexistence of the unitary evolution of the ainphtude to- 
gether with the asymptotic value of the GCD is a striking 
result about the behavior of the system. 

We study the entanglement between the coin and the 
position in the QW on the line and we show that the 
behavior of the entropy of entanglement depends on the 
GCD. We also show that the asymptotic entanglement 
is maximized when the evolution of the GCD follows a 
Markovian process. However the entropy of entangle- 
ment does not describe correctly the degree of disorder 
of the GCD; this is well described by the Shannon en- 
tropy. In previous works [TTI - fTsl [2l| the dependence 
of the asymptotic entropy of entanglement with the ini- 
tial conditions was studied, here we provide an analyti- 



cal recipe to obtain a predetermined entanglement using 
extended Gaussian initial conditions. In other words, 
starting from a given value of the entropy of entangle- 
ment it is possible to choose the corresponding initial 
conditions. These exact expressions can be also used to 
obtain a predetermined asymptotic GCD, that is, start- 
ing from a given asymptotic limit of the GCD obtain the 
corresponding initial conditions for the QW. 
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and Carlos Navarrete Benlloch and the support from 
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